CPT Violating Decoherence and LSND: a possible window to Planck scale Physics 
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Decoherence has the potential to explain all existing neutrino data including LSND results, with- 
out enlarging the neutrino sector. This particular form of CPT violation can preserve the equality of 
masses and mixing angles between particle and antiparticle sectors, and still provide seizable differ- 
ences in the oscillation patterns. A simplified minimal model of decoherence can explain the existing 
neutrino data as well as the standard three oscillation scenario, while making dramatic predictions 
for the upcoming experiments. Such a model can easily accomodate the LSND result but cannot 
fit the spectral distortions seen by KamLAND. Some comments on the order of the decoherence 
parameters in connection with theoretically expected values from some models of quantum-gravity 
are given. In particular, the quantum gravity decoherence as a primary origin of the neutrino mass 
differences scenario is explored, and even a speculative link between the neutrino mass-difference 
scale and the dark energy density component of the Universe today is drawn. 



INTRODUCTION: THE LSND PUZZLE AND 
ABANDONING OF SYMMETRIES 



Although the Special Theory of Relativity, a theory of 
fiat space time physics based on Lorentz symmetry, is 
very well tested, and in fact next year it celebrates a cen- 
tury of enormous success, having passed very stringent 
experimental precision tests, a quantum theory of Grav- 
ity, that is a consistent quantized version of Einstein's 
General Relativity, a dynamical model for curved space 
times, still eludes us. The main reason for this is the 
lack of any concrete observational evidence on the struc- 
ture of space time at the characteristic scale of quantum 
gravity (QG), the (four-dimensional) Planck mass scale 
Mp = 10 19 GeV. The situation concerning QG is to be 
contrasted with that characterizing its classical counter- 
part, General Relativity, and the quantum version of Spe- 
cial Relativistic field theories, both of which can be tested 
to enormous precision to date. For instance, the perihe- 
lion precession of planets such as mercury, as well as the 
deflection of light by the sun, made the General Theory 
of Relativity an instant success, almost immediately after 
it was proposed in 1915. Similarly, the modern versions 
of relativistic quantum field theories, including the foun- 
dations of the standard model of particle physics, which 
have been tested extremely well up to now, leave little 
doubt on the validity of Special Relativity as a quantum 
theory of flat space time, that is in the absence of gravi- 
tational effects. 

One would hope that, like any other successful physi- 
cal theory, a physically relevant quantum theory of grav- 
ity should lead to experimental predictions that should 
be testable in the foreseeable future. In the QG case, 
however, such predictions may be subtle, due to the 
extremely weak nature of the gravitational interaction 
as compared with the rest of the known interactions 
in nature. Indeed, any hope for a true "phenomenol- 



ogy of QG" may be easily proven to be wishful think- 
ing, if all the symmetries and properties that character- 
ize the current version of what is called particle physics 
phenomenology are still valid in a full quantum the- 
ory of gravity. This is mainly due to the fact that the 
dimension-full coupling constant of gravity, the Newton 
constant Gn = l/Mp, appears as a very strong sup- 
pression factor of any physical observable that could be 
associated with predictions of quantum gravity, in a the- 
ory in which Lorentz invariance and the laws of quantum 
mechanics (or better quantum field theory), such as uni- 
tary temporal evolution and locality of interactions, still 
hold. 

In recent years, however, more and more physicists 
contemplate the idea that such laws, which are charac- 
teristic of a flat space time quantum field theory, may 
not be valid in a full theory of curved space times. For 
instance, the CPT theorem, one of the most profound re- 
sults of quantum field theory , which is a consequence 
of Lorentz invariance, locality, as well as quantum me- 
chanics (specifically unitary evolution of a system), may 
not characterize a quantum gravity theory. The possibil- 
ity of a violation of CPT invariance by quantum gravity 
has been raised in a number of theoretical models of QG 
that go beyond conventional local quantum field theoretic 
treatments of gravity 0, 0, IE @ • 

In a phenomenological, rather model-independent set- 
ting, CPT Violation has been invoked recently in 
an attempt to explain within a three generation sce- 
nario, without the introduction of sterile neutrinos, a 
puzzling experimental result in neutrino physics, namely 
the claims from LSND Collaboration 7] on evidence of 
oscillations in the antineutrino sector —* v e , due to 
observed excess of v e events, but lack of evidence (de- 
pending though on interpretation) for oscillations in the 
corresponding neutrino sector — > v e . In order to ex- 
plain these results within the conventional oscillation see- 
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nario, one should invoke different neutrino mass differ- 
ences in the particle and antiparticle sectors of the model. 
This would signify CPT violation. Specifically, from so- 
lar and atmospheric neutrino data, the mass squared dif- 
ferences that account for the observed oscillations are of 
order Am 2 ~ 10~ 5 eV 2 and 10~ 3 eV 2 , respectively, while 
the LSND result would require Am 2 ~ 10 _1 eV 2 , which 
would explain the suppressed signal in the neutrino sec- 
tor and the strong signal in the antineutrino one. To 
reconcile, therefore, the LSND result with the rest of 
the available neutrino data, in conventional oscillation 
scenaria, one would have to invoke CPT violating mass 
differences between neutrinos and antineutrinos of order 
|m| - m 2 | ~ 1CT 1 eV 2 . Although such a CPT viola- 
tion may look rather drastic, the most stringent limit 
we have on the CPT symmetry today, the one com- 
ing from the neutral kaon system, if reinterpreted as a 
limit on the possible differences in mass squared, gives: 
\m 2 {K°) -m 2 (K°)\ < 0.25 eV 2 showing that such differ- 
ences are just barely explored even in the kaon system. 

Subsequent global analyses of neutrino data , how- 
ever, seem to disfavor two and three generation CPT vi- 
olating scenaria 0, thereby leaving the CPT violat- 
ing scenaria involving sterile neutrinos as the only sur- 
viving possibility for the explanation of the LSND re- 
sults. Specifically, it has been argued in [Toj| . that al- 
lowing for CPT violation, the mixing matrix elements 
with the forth generation between neutrino and antineu- 
trino sectors need no longer be the same, thereby evad- 
ing stringent experimental constraints that killed three- 
generation models. 

Is, however, CPT violation evidenced only as an in- 
equality between neutrino-antineutrino masses the only 
way a violation of this symmetry can manifest itself in 
nature? Such a question becomes extremely relevant for 
the case of LSND, because it is possible that other mech- 
anisms leading to CPT violation exist, unrelated, in prin- 
ciple, to mass differences between particles and antiparti- 
cles. Such additional mechanisms for CPT violation may 
well be capable of explaining the LSND results within a 
three generation scenario without invoking a sterile neu- 
trino (a scenario which, on the other hand, is getting 
totally excluded as new experimental data become avail- 
able). It is therefore necessary to explore whether alter- 
native ways exist to account for the LSND result without 
invoking extra (sterile) neutrino states. 

As we shall argue in this article, quantum decoherence 
may be the key to answer this question. Indeed, quantum 
decoherence in matter propagation occurs when the mat- 
ter subsystem interacts with an 'environment', according 
to the rules of open-system quantum mechanics. At a 
fundamental level, such a decoherence may be the result 
of propagation of matter in quantum gravity space-time 
backgrounds with 'fuzzy' properties, which may be re- 
sponsible for violation of CPT in a way not necessarily 
related to mass differences between particles and antipar- 



ticles. 

It is the point of this work to claim that decoherence 
scenaria combined with conventional oscillations, led as 
usual by the mass differences between the mass eigen- 
states, can account for all the available neutrino data, 
including the LSND result, within a three generation 
model. The LSND result would then evidence CPT vi- 
olation in the sense of different decohering interactions 
between particle and antiparticle sectors, while the mass 
differences between the two sectors remain the same. For 
detailed reviews of decoherence in two-level systems, with 
emphasis on phenomenology including neutrinos, we re- 
fer the reader to the literature I n what follows we 
shall only expose the very basic features, which we shall 
make use of in this work in order to perform our phe- 
nomenological analysis. It is important to stress that in 
decoherence scenaria CPT symmetry is violated in its 
strong form, in the sense of the CPT operator not being 
well defined. This is an important issue which we now 
come to discuss briefly. 

QUANTUM GRAVITY AND DECOHERENCE: 
BRIEF REVIEW 

A characteristic example of such a violation occurs in 
quantum gravity models that involve singular space-time 
configurations, integrated over in a path integral formal- 
ism, which are such that the axioms of quantum field the- 
ory, as well as conventional quantum mechanical behav- 
ior, cannot be maintained 3] . Such configurations consist 
of wormholes, microscopic (Planck size) black holes, and 
other topologically non-trivial solitonic objects, such as 
geons 12], etc. Collectively, we may call such config- 
urations space time foam, a terminology given by J. A. 
Wheeler [l3| who first conceived the idea that the struc- 
ture of the quantum space time at Planck scales (10~ 35 
m) may indeed be fuzzy. 

It has been argued that, as result, a mixed state de- 
scription must be used (QG-induced decoherence) Q-El- 
given that such objects cannot be accessible to low- 
energy observers, and as such must be traced over in 
an effective field theory context. For the case of micro- 
scopic black holes this can be readily understood by the 
loss of information across microscopic event horizons, as 
a corollary of which it has been argued Q| > already back 
in 1978, that CPT invariance in its strong form must be 
abandoned in a foamy quantum gravity theory. Such a 
breakdown of CPT symmetry is a fundamental one, and, 
in particular, implies that a proper CPT operator may 
be ill defined in such QG decoherence cases. 

The main reason behind such a failure is the non fac- 
torisability of the super-scattering matrix, $j_connecting 
asymptotic "in" and "out" density matrices |3j. The lat- 
ter describe mixed quantum states due to the existence 
of a quantum-gravitational "environment" of microscopic 



3 



singular space-time configurations. The presence of the 
environment leads to decoherence, like in any other case 
of an open quantum-mechanical system, implying an evo- 
lution from an initial pure quantum mechanical state 
|0 > to a mixed state p^ = Tr|-0 >< (the trace is 
over unobserved degrees of freedom, 'lost' in the environ- 
ment). This in turn implies 14] the non invertibility of 
the $-matrix, as a quantum mechanical operator acting 
on density matrices. The associated CPT non invariance 
is then linked to the fact that if CPT were a well defined 
operator, stemming from the invariance of the system 
under its action, then $ would have been invertible, in 
obvious contradiction with the information loss and de- 
coherence inherent to the problem. 

As emphasized recently in 0] this would also imply 
that the concept of antiparticle, may be ill defined, with 
interesting testable consequences in multiparticle situa- 
tions, such as those encountered in (EPR-like entangled) 
meson factories. Such CPT violation is distinct from 
other violations of this symmetry as a result of Lorentz 
symmetry breakdown in the so-called Standard Model 
extension (SME) 0, or locality violations which have 
been considered in the recent literature as other exper- 
imentally testable approaches to quantum gravity. It 
should be noted that these latter violations of symmetries 
may or may not coexist with QG decoherence, given that 
Lorentz invariance is not necessarily inconsistent with de- 
coherence 

Some caution should be paid regarding CPT Viola- 
tion through decoherence. As emphasized in [J^j, from a 
formal view point, the non- invertibility of the $-matrix, 
which implies a strong violation of CPT in the sense de- 
scribed above, does not preclude a softer form of CPT 
invariance, in the sense that any strong form of CPT vio- 
lation does not necessarily have to show up in any single 
experimental measurement. This implies that, despite 
the general evolution of pure to mixed states, it may still 
be possible in the laboratory to ensure that the system 
evolves from an initial pure state ip to a single final state 
<j), and that the weak form of CPT invariance is mani- 
fested through the equality of probabilities between the 
states -0, 4>: 

p^^cp) = p{e- 1 cj ) ^e^) , (i) 

with 9 a CPT operator acting on the subspace of pure 
state vectors in the laboratory, such that for density ma- 
trices Op = 0p0\ (9+ = -9-\ and $t = e^G" 1 , but 
$ ' ^$ _1 J unless full CPT invariance holds, which is 
the case of no decoherence. If this is the case, then the 
decoherence-induced CPT violation will not show up in 
any experimental measurement. It is therefore impor- 
tant to check whether it is possible to test experimentally 
which case is realized in nature. Fortunately this can be 
done in a clear way. 

Indeed, the main consequence of the quantum-gravity 
induced decoherence is the modified non-quantum me- 



chanical evolution of a matter system under considera- 
tion, as a result of the interaction with the QG environ- 
ment 0. Such an evolution will result in many observ- 
able consequences, which may be testable in the near fu- 
ture with a precision that in some cases can reach Planck 
scale sensitivity. Historically, the first probe of such non- 
quantum mechanical decoherence and CPT violatin g be - 
havior (in the sense of the strong form violation of 
is the neutral kaon system 0, 0| , in both single-particle 
experiments, such as the CPLEAR 0] and KTeV exper- 
iments as well as multi-particle situations, such as 
kaon ((/>) factories [l^. Other neutral mesons, such as B 
mesons [i^, can also be used as sensitive probes of QG 
decoherence. 

In single-particle neutral meson experiments one looks 
at time profiles of asymmetries 0. The mass difference 
between the energy eigenstates Kl, Ks is responsible for 
the 'bulk' behavior of the asymmetry, for instance lead- 
ing to an oscillatory pattern in some of the asymmetries, 
such as A^, decoherence on the other hand can be re- 
sponsible for a slight distortion of such a behavior 0]. 
The non observation of such distortions results in strin- 
gent bounds on the decoherence parameters. Mathemat- 
ically, one does not have to have a detailed knowledge 
of the dynamics of quantum gravity in order to perform 
experimental low-energy tests of decoherence. This can 
be achieved following the so-called Lindblad or mathe- 
matical semi-groups approach to decoherence |2l| , which 
is a very efficient way of studying open systems in quan- 
tum mechanics. The time irreversibility in the evolution 
of such semigroups, which is linked to decoherence, is 
inherent in the mathematical property of the lack of an 
inverse in the semigroup. This approach has been fol- 
lowed for the study of quantum-gravity decoherence in 
the case of neutral kaons in |i|4|. 

In the parameterization of [2j for the decoherence ef- 
fects, using three decoherence parameters with dimen- 
sions of energy, a, j3, 7, one finds the following bounds 
from the CPLEAR experiment [T^: 

a < 4 x 1(T 17 GcV, \j3\ < 2.3 x 1(T 19 GeV, 

7 < 3.7 x 1(T 21 GeV (2) 

where the positivity of p, required by the fact that its 
diagonal elements express probability densities, implies 
a, 7 > 0, and > /3 2 . 

The Lindblad approach to decoherence does not re- 
quire any detailed knowledge of the environment, apart 
from energy conservation, entropy increase and complete 
positivity of the (reduced) density matrix p{t) of the sub- 
system under consideration. The basic evolution equa- 
tion for the (reduced) density matrix of the subsystem in 
the Lindblad approach is linear in p(t) and reads: 

3 

(3) 
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where H c g is the effective Hamiltonian of the subsystem, 
and the operators bj represent the interaction with the 
environment, and are assumed bounded. Notice that the 
Lindblad part cannot be written as a commutator (of a 
Hamiltonian function) with p. Environmental contribu- 
tions that can be cast in Hamiltonian evolution (commu- 
tator form) are absorbed in H c q. 

It must be noted at this stage that the requirement 
of complete positivity, which essentially pertains to the 
positivity of the map p(t) as the time evolves in the case 
of many particle situations, such as meson factories (two- 
kaon states (0- factory), or two- B- meson states etc.), may 
not be an exact property of quantum gravity, whose in- 
teractions with the environment could be non linear |22| . 
Nevertheless, complete positivity leads to a convenient 
and simple parametrization, and it has been assumed 
so far in many phenomenological analyses of quantum 
gravity decoherence in generic two state systems, such 
as two- flavor neutrino systems In fact, in the 

parametrization of Q, the imposition of complete posi- 
tivity leads to a = /3 = 0, leaving only 7 > as the only 
decoherence parameter in the two-state system |2Fi| . 

Formally, the bounded Lindblad operators of an TV- 
level quantum mechanical system can be expanded in a 
basis of matrices satisfying standard commutation rela- 
tions of Lie groups. For a two-level system 0, such 
matrices are the SU(2) generators (Pauli matrices) plus 
the 2x2 identity operator, while for a three level sys- 
tem |2(il |. which will be relevant for our purposes in this 
article, the basis comprises of the eight Gell-Mann SU (3) 
matrices Aj , £ = 1, ... 8 plus the 3x3 identity matrix 

13x3- 

Let Jn, fi = 0, . . . 8(3) be a set of SU(3) (SU(2)) gener- 
ators for a three (two)-level system; then, one may expand 
the various terms in J^J in terms of to arrive at the 
generic form: 

ij V 

/i,^ = 0,...iV 2 -l, i,j = l,...N 2 -l (4) 

with N = 3(2) for three(two) level systems, and fijk 
the structure constants of the SU(N) group. The re- 
quirement for entropy increase implies the hermiticity 
of the Lindblad operators as well as the fact that 
the matrix C of the the non-Hamiltonian part of the 
evolution has the properties that Cq^ = C^q = 0, 
Aj = ISfc/.m^&fc ] hmkfik], with the notation bj = 

In the two-level case of |2j the decoherence matrix £ M „ 
is parametrized by a 4 x 4 matrix, whose non vanishing 
entries are occupied by the three parameters with the di- 
mensions of energy a, (3, 7 with the properties mentioned 
above. If the requirement of a completely positive map 
p(t) is imposed, then the 4x4 matrix C becomes diag- 



onal, with only one non vanis hing entry occupied by the 
decoherence parameter 7 > [25j. 

TWO GENERATION NEUTRINO MODELS AND 
DECOHERENCE: BRIEF REVIEW 

If the above formalism is applied to two-level neutrino 
oscillation physics, in an attempt to discuss the influence 
of a quantum-gravity environment on neutrino oscilla- 
tions, then there are two possible forms of the matrix C: 
C = Diag (0, 0, —7, 0) in the case where energy and lepton 
number are conserved, and C = Diag (0, 0, 0, —7) in the 
case that energy and lepton number are violated by quan- 
tum gravity, but flavor is conserved. In the energy 
conserving case has been considered, with the conclusion 
that pure decoherence, i.e. induced oscillations in the ab- 
sence of neutrino mass differences, was incompatible with 
the observed phenomenology of neutrino oscillations. 

This is in agreement with the kaon case , where again 
the observed bulk behavior of the time profiles of the var- 
ious kaon asymmetries, as well as the observed CP viola- 
tion, cannot be explained solely by quantum decoherence, 
the latter providing only distortions to the 'bulk' behav- 
ior. Including the mass-induced oscillations, combined 
with decoherence, and comparing with atmospheric neu- 
trino data, the authors of [23| managed to obtain bounds 
for the decoherence parameter 7 in the presence of mass 
differences (oscillations); as in the kaon case Q, such 
bounds are obtained merely by the non-observation of 
the slight distortions to the mass-induced neutrino oscil- 
lations that the presence of decoherence would cause. 

The so-obtained limits for the decoherence Lindblad 
parameter 7 in can be summarized as follows: adopt- 
ing the parametrization 7 = 7 (-E/GeV) n , they consid- 
ered three cases: (a) n=0, in which the experimental data 
imply 70 < 10 -23 GeV, (b) n=2, which is the case ex- 
pected in some rather optimistic (from the prospect of 
experimental detection) theoretical models of quantum 
gravity 0, , for which the most stringent bound de- 
rived is 70 < 0.9 x 1CP 27 , and (c) n=-l, which is the 
case in which the decoherence parameter exhibits an en- 
ergy dependence that mimics the conventional oscillation 
scenario, for which 70 < 2 x 1CP 21 GeV 2 . 

It is worth noticing, that as far as case (b) is concerned, 
at least by dimensional considerations, one is tempted to 
assume theoretical values of the decoherence parameter 
which are of the form E 2 /Mqg, with Mqq the effec- 
tive quantum gravity scale 'felt' by the neutrinos. The 
above bounds indicate a value Mqq ~ 10 27 GeV for 
E = C(few GeV) which are typical in atmospheric neu- 
trino experiments. This is much higher than the theo- 
retically expected (on naturalness grounds) Planck value 
Mqq ~ 10 19 GeV. This is to be contrasted with the case 
of kaons Q, where, as mentioned above, the bound on 
the decoherence parameter 7 < 10 _21 GeV obtained from 
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CPLEAR [lTj , was found to be much closer to the Planck 
scale, the natural QG scale. On the other hand, the case 
(c), with an inverse energy E dependence, is purely phe- 
nomenological and does not stem directly from a specific 
model of decoherence. The bound on 70 < 10~ 23 (GeV) 2 
obtained in that case is of the order of the quantity 
Am 2 /E for atmospheric neutrino energies, and square 
mass differences of order 10 -5 eV 2 which is assumed to 
be the right order for solar neutrino oscillations in the 
conventional oscillation scenario. 

Unfortunately, the analysis of the data presented in 
|23| within a two-generation scenario did not show any 
concrete evidence for such a type-(c) decoherence, even 
after inclusion [2?J of the first K2K data Eii. in com- 
bination with Super Kamiokande (SK) data .29]. This 
prompted these authors to draw the pessimistic conclu- 
sion that the extension of the neutrino oscillation plus 
decoherence scenario to the fully-fledged case of three 
neutrino generations was not worthy. Adding to this pes- 
simism, there were also some theoretical estimates of the 
decoherence parameter in specific models of quantum- 
gravi ty i nduced decoherence for generic two-level sys- 
tems (33, according to which 70 ~ (Am 2 ) 2 /E 2 Mp, with 
Am 2 the mass-squared difference between, say, the two 
neutrinos (or kaons etc). For Mp ~ 10 19 GeV, the Planck 
scale, this is beyond any prospects for experimental de- 
tection in oscillation experiments in the foreseeable fu- 
ture. However, as remarked in 122 1, s uch estimates are 
specific to the model considered in |30j , and indeed there 
are concrete examples of theoretical (stringy) models of 
quantum-gravity-induced decoherence where they are not 
applicable. We shall discuss such issues later on in our 
article. 



EXTENSION TO THREE GENERATIONS: 
COMBINING OSCILLATIONS WITH 
DECOHERENCE 

Pessimistic thoughts have also appeared in the work 
of [2(|, which considered the extension of the com- 
pletely positive decoherence scenario to the standard 
three-generation neutrino oscillations case. The above- 
described three-state Lindblad problem has been adopted 
for the discussion of this case. The relativistic neutrino 
Hamiltonian H c s ~ p 2 + m 2 /2p, with m the neutrino 
mass, has been used as the Hamiltonian of the subsys- 
tem in the evolution of eq. © . 

In terms of the generators J^, fj, = 0,...8 of 
the SU(3) group, H c g can be expanded as [26| : 

H cS = l p ^W (6P 2 + Eli mfj Jo + ^(Am 2 2 )J3 + 
2 ^ (Am 2 3 + Am| 3 ) with the obvious notation 
Am 2 ,- =m 2 -m 2 , i,j = 1,2,3. 

The analysis of [26| assumed ad hoc a diagonal 
form for the 9x9 decoherence matrix C in 



[£„„] = Diag (0, -71, -72, -73, -74, -75, -76, -77, ~7s) 
in direct analogy with the two-level case of complete pos- 
itivity [23|,|25|]. As we have mentioned already, there is no 
strong physical motivation behind such restricted forms 
of decoherence. This assumption, however, leads to the 
simplest possible decoherence models, and, for our phe- 
nomenological purposes in this work, we will assume the 
above form, which we will use to fit all the available neu- 
trino data. It must be clear to the reader though, that 
such a simplification, if proven to be successful (which, as 
we shall argue below, is the case here), just adds more in 
favor of decoherence models, given the restricted num- 
ber of available parameters for the fit in this case. In 
fact, any other non-minimal scenario will have it easier 
to accommodate data because it will have more degrees 
of freedom available for such a purpose. 

Specifically we shall look at transition probabili- 
ties 0]: 

i + i^e^^-i^^f ( 5 ) 

where a, j3 = e, fj,,T stand for the three neutrino flavors, 
and Latin indices run over 1, ... 8. The quantities are 
the eigenvalues of the matrix M. appearing in the evolu- 
tion ijlf. after taking into account probability conserva- 
tion, which decouples po (t) = \/2/3, leaving the remain- 
ing equations in the form: dpk/dt = ^ZjMkjPj- The 
matrices X>y are the matrices that diagonalize M. |2lj . 
Explicit forms of these matrices, the eigenvalues Afe, and 
consequently the transition probabilities © , are given in 

0. 

The important point to stress is that, in generic models 
of oscillation plus decoherence, the eigenvalues Afe depend 
on both the decoherence parameters 7; and the mass 
differences Amf,-. For instance, Ai = 5 [—(71 +72) — 
\J (72 — 71 ) 2 — 4A 2 2 ], with the notation Ay- = Am 2 J /2p, 
i, j = 1,2,3. Note that, to leading order in the (small) 
squared-mass differences, one may replace p by the total 
neutrino energy E, and this will be understood in what 
follows. For detailed expressions of the rest of the pa- 
rameters we refer the reader to |2f|. However, we note, 
for future use, that it is a generic feature of the Afe to 
depend on the quantities fiy which are given by 



fil2 = 


\/(72 


-7i) 2 


-4A? 2 


fil3 = 


\/(75 


-74) 2 


-4A 2 3 


^23 = 


J{l7 


-7e) 2 


-4A 2 3 



From the above expressions for the eigenvalues Afe, it be- 
comes clear that, when decoherence and oscillations are 
present simultaneously, one should distinguish two cases, 
according to the relative magnitudes of and A7fe; = 
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| > |A 7fc£ |, and (ii) 2|Ay| < |A 7w |. In 
the former case, the probabilities JSJ) contain trigonomet- 
ric (sine and cosine) functions, whilst in the latter they 
exhibit hyperbolic sin and cosine dependence. 

Assuming mixing between the flavors, amounts to ex- 
pressing neutrino flavor eigenstates \v a >, a = e, p,,T in 
terms of mass eigenstates \Vj, >, i = 1,2,3 through a 
(unitary) matrix U: \v a >= X)i=i Uai\ v i - > - This im- 
plies that the density matrix of a flavor state p a can be 
expressed in terms of mass eigenstates as: p a = \v a >< 
u a \ = J2i j UaiUajWi >< Vj\- From this we can deter- 
mine = 2Tr(p Q l 7 AI ), a quantity needed to calculate the 
transition probabilities |J5J. 

The important comment we would like to raise at this 
point is that, when considering the above probabilities in 
the antineutrino sector, the respective decoherence pa- 
rameters 7i in general may be different from the corre- 
sponding ones in the neutrino sector, as a result of the 
strong form of CPT violation. In fact, as we shall discuss 
next, this will be crucial for accommodating the LSND 
result without conflicting with the rest of the available 
neutrino data. This feature is totally unrelated to mass 
differences between flavors. 



COMBINING DECOHERENCE, OSCILLATIONS 
AND THE LSND RESULT 



In |26j a pessimistic conclusion was drawn on the "clear 
incompatibility between neutrino data and theoretical 
expectations", as followed by their qualitative tests for 
decoherence. It is a key feature of our present work to 
point out that we do not share at all this point of view. 
In fact, as we shall demonstrate below, if one takes into 
account all the available neutrino data, including the fi- 
nal LSND results Q, which the authors of [2f| did not 
do, and allows for the above mentioned CPT violation in 
the decoherence sector, then one will arrive at exactly the 
opposite conclusion, namely that three-generation deco- 
herence and oscillations can fit the data successfully! 

As we shall argue here, compatibility of all avail- 
able data, including CHOOZ and LSND, can be 
achieved through a set of decoherence parameters jj 
with energy dependences "JjE and jj/E, with 7? ~ 
10~ 18 ,10~ 24 (GeV) 2 , respectively, for some j's. As it 
will be clear later, the quality of the fit diminishes to the 
level of the standard three-oscillation case once Kam- 
LAND spectral distortion data js^l are included but is 
significanlty higher otherwise, showing that either addi- 
tional energy dependences are needed or non-diagonal 
elements are required to get better agreement with data. 

Some important remarks are in order. First of all, 
in the analysis of |2y] pure decoherence is excluded in 
three-generation scenaria, as in two generation ones, due 
to the fact that the transition probabilities in the case 
Am 2 - = (pure decoherence) are such that the survival 



probabilities in both sectors are equal, i.e. P(y n —> v a ) 
P{y a ^v a ). From (J5J we have in this case [26(: 
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From the CHOOZ experiment [31] . for which L/E ~ 
10 3 /3 m/MeV, we have that (Pv e ^v e ) - 1, while the 
K2K experiment [H with L/E ~ 250/1.3 km/GcV has 
observed events compatible with (P Vfl _n, ) — 0.7, thereby 
contradicting the theoretical predictions J7J) of pure de- 
coherence. 

However, this conclusion is based on the fact that in 
the antineutrino sector the decoherence matrix is the 
same as that in the neutrino sector. In general this need 
not be the case, in view of CPT violation, which could 
imply a different interaction of the antiparticle with the 
gravitational environment as compared with the particle. 

In our tests we take into account this possibility, but 
as can be seen from the figure, and will be evident later 
on, once the analysis will be explained in some depth, 
pure decoherence can be excluded also in this case, as it 
is clearly incompatible with the totality of the available 
data. 

In order to check our model, we have performed a x 2 
comparison (as opposed to a \ 2 fit) to SuperKamiokande 
sub-GeV and multi GeV data (the forty data points 
that are shown in the figures), CHOOZ data (15 data 
points) LSND (1 datum) and KamLAND spectral infor- 
mation |32l | (sampled in 13 bins), for a sample point in 
the vast parameter space of our extremely simplified ver- 
sion of decoherence models. Let us emphasize that we 
have not performed a x 2_n t and therefore the point we 
are selecting (by "eye" and not by \) is n °t optimized 
to give the best fit to the existing data. Instead, it must 
be regarded as one among the many equally good sons in 
this family of solutions, being extremely possible to find 
a better fitting one through a complete (and highly time 
consuming) scan over the whole parameter space. 

Cutting the long story short, and to make the analysis 
easier, we have set all the 7$ in the neutrino sector to 
zero, restricting this way, all the decoherence effects to 
the antineutrino one where, we have assumed for the sake 
of simplicity, 

7s = li+i f° r » = 1, 4, 6, 7 and 71 = 74 , 73 = 78 (8) 

Furthermore, we have also set the CP violating phase of 
the NMS matrix to zero, so that all the mixing matrix 
elements become real. 

With these assumptions, the otherwise cumbersome 
expression for the transition probability for the antineu- 
trino sector takes the form, 
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where the Ojj were defined in the previous section and 
are the same in both sectors (due to our choice of "fi's) 
and 



p^ = 2Re(U* 1 U a2 ) 
P^ = \U al \ 2 - \U a2 \ 2 
P 2 = 2Re(U* al U a3 ) 
p% = 2 Re(U* 2 U a3 ) 

Ps 



(10) 
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where the mixing matrices are the same as in the neutrino 
sector. For the neutrino sector, as there are no decoher- 
ence effects, the standard expression for the transition 
probability is valid. 

It is obvious now that, since the neutrino sector does 
not suffer from decoherence, there is no need to include 
the solar data into the fit. We are guaranteed to have an 
excellent agreement with solar data, as long as we keep 
the relevant mass difference and mixing angle within the 
LMA region, something which we shall certainly do. 

As mentioned previously, CPT violation is driven by, 
and restricted to, the decoherence parameters, and hence 
masses and mixing angles are the same in both sectors, 
and selected to be 

= ArrTH 2 = 7 • 10~ 5 eV 2 , 
^_Arnii 2 = 2.5 • KT 3 eV 2 , 

#23 = 6>23 = TT/4^12 = 012 = .45, 
013 = 013 = -05, 

as indicated by the state of the art analysis. 

For the decoherence parameters we have chosen (c.f. 



Am 2 2 
Am^ 3 = 



)) 



and 



73 = 76 = 77 



7^=i.io- 2 7£, 



(11) 



where E is the neutrino energy, and barred quantities 
refer to the antineutrinos, given that decoherence takes 
place only in this sector in our model. All the other 
parameters are assumed to be zero. All in all, we have 
introduced only two new parameters, two new degrees 
of freedom, 71 and 73, and we shall try to explain with 
them all the available experimental data. 

It can be checked straightforwardly, by means of 
©, that for the regime of the parameters relevant to 
the experiments we are considering in this work, this 
parametrization guarantees the positivity of the relevant 
probabilities. This is an important issue, because usu- 
ally negative probabilities are viewed as a signal of in- 
consistency in parameterizing the pertinent decoherence 



effects. In our case, with two different energy depen- 
dences in the (non constant) decoherence coefficients, 
there seems to be no master condition that guarantees 
the complete positivity in a way independent of the ener- 
gies/momenta of the neutrinos. As we shall discuss later 
on in the article, we attribute the l/-E-dependent deco- 
herent coefficients to conventional matter effects, while 
the ^-dependent ones are associated with novel effects of 
quantum gravity, which increase with the energy of the 
probe. The negative probabilities that may occur out- 
side the regime of parameters to be specified below are 
then interpreted as implying simply that our linear (and 
simplified) parametrization of quantum gravity decoher- 
ence QUI) ceases to be valid in such regimes, and more 
complicated, probably non linear, entanglement may be 
in operation, as expected in a quantum theory of grav- 
ity 0. 

With these in mind, we now observe that, in our simpli- 
fied model of decoherence, parametrized by a diagonal de- 
coherence matrix, positivity of the relevant probabilities 
seems to be guaranteed for ^f x L > 7 3 L, which, with our 
parametrization JH), implies: L/E < 10 24 -10 25 GcV~ 2 , 
and E > 1 MeV. These sufficient conditions are met by 
all the current and planned terrestrial (anti)neutrino ex- 
periments. Outside this regime, our parametrization sim- 
ply fails, and one needs to resort to more complicated sit- 
uations, which fall beyond the scope of the present paper, 
and will be presented elsewhere. 

At this point it is important to stress that the inclu- 
sion of two new degrees of freedom is not sufficient to 
guarantee that one will indeed be able to account for all 
the experimental observations. We have to keep in mind 
that, in no-decoherence situations, the addition of a ster- 
ile neutrino (which comes along with four new degrees of 
freedom -excluding again the possibility of CP violating 
phases) did not seem to be sufficient for matching all the 
available experimental data, at least in CPT conserving 
situations. 

In order to test our model with these two decoher- 
ence parameters in the antineutrino sector, we have cal- 
culated the zenith angle dependence of the ratio "ob- 
served/ (expected in the no oscillation case)", for muon 
and electron atmospheric neutrinos, for the sub-GeV and 
multi-GeV energy ranges, when mixing is taken into ac- 
count. Since matter effects are important for atmospheric 
neutrinos, we have implemented them through a two- 
shell model, where the density in the mantle (core) is 
taken to be roughly 3.35 (8.44) gr/cm 3 , and the core ra- 
dius is taken to be 2887 km. We should note at this stage 
that a "fake" CPT Violation appears due to matter ef- 
fects, arising from a relative sign difference of the matter 
potential between the respective interactions of neutrinos 
and antineutrinos with ordinary matter. This, however, 
is easily disentangled from our genuine (due to quantum 
gravity) CPT Violation, used here to parametrize our 
model fit to LSND results; indeed, a systematic study of 
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such effects 33] has shown that "fake" CPT Violation in- 
creases with the oscillation length, but decreases with the 
neutrino energy, E, vanishing in the limit E — > oo; more- 
over, no independent information regarding such effects 
can be obtained by looking at the antineutrino sector, as 
compared with data from the neutrino sector, due to the 
fact that in the presence of "fake" CPT Violation, but 
in the absence of any genuine CPT breaking, the perti- 
nent CPT probability differences between neutrinos and 
antincutrinos arc related 



ApCPT 



= -AP§I T , where 

pa ' 



P-5 and the Greek indices denote neu- 



a/3 * p^l 

trino flavors. These features are to be contrasted with 
our dominant decoherence effects 7j (flip . proportional 
to the antineutrino energy, E, which are dominant only 
in the antineutrino sector. For the same reason, our ef- 
fects can be disentangled from "fake" decoherence effects 
arising from Gaussian averages of the oscillation proba- 
bility due to, say, uncertainties in the energy of the neu- 
trino beams |34| . which are the same for both neutrinos 
and antineutrinos. We, therefore, claim that the com- 
plex energy dependence in , with both L ■ E and L / E 
terms being present in the antineutrino sector, may be a 
characteristic feature of new physics, with the L ■ E terms 
being related to quantum-gravity induced (genuine) CPT 
Violating decoherence. 

The results are shown in Fig. 1 (c), where, for the 
sake of comparison, we have also included the experi- 
mental data. We also present in that figure the pure 
decoherence scenario in the antineutrino sector (a), as 
well as in both sectors (b). For completeness, we also 
present a scenario with neutrino mixing but with deco- 
herence operative in both sectors (d). The conclusion is 
straightforward: while pure decoherence appears to be 
excluded, decoherence plus mixing provides an astonish- 
ing agreement with experiment |45| . 

As bare eye comparisons can be misleading, we have 
also calculated the \ 2 value for each of the cases, defining 
the atmospheric % 2 as 
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Here a a ^ are the statistical errors, the ratios R a ,i be- 
tween the observed and predicted signal can be written 
as 



R 



exp 



MC 



(13) 



(with a indicating the lepton flavor and i counting the 
different bins, ten in total) and M, S stand for the multi- 
GeV and sub-GeV data respectively. For the CHOOZ 
experiment we used the 15 data points with their statis- 
tical errors, where in each bin we averaged the probability 
over energy and for LSND one datum has been included. 
The results with which we hope all our claims become 
crystal clear are summarized in Table 1 and 2 , were we 





FIG. 1: Decoherence fits, from top to bottom: (a) pure deco- 
herence in antineutrino sector, (b) pure decoherence in both 
sectors, (c) mixing plus decoherence in the antineutrino sector 
only, (d) mixing plus decoherence in both sectors. The dots 
correspond to SK data. 



present the \ 2 comparison for the following cases: (a) 
pure decoherence in the antineutrino sector, (b) pure de- 
coherence in both sectors, (c) mixing plus decoherence in 
the antineutrino sector, (d) mixing plus decoherence in 
both sectors, and (e) mixing only - the standard scenario 
for all experimental data with (Table 2) and without (Ta- 
ble 1) KamLAND spectral distortion . 

From the tables it becomes clear that the mixing plus 
decoherence scenario in the antineutrino sector can easily 
account for all the available experimental information, in- 
cluding LSND when KamLAND spectral distortion data 
are excluded but does not significantly improve over the 
standard case once they are included, showing clearly the 
limitations of our simplified scenario. 

Such a change can be easily undestood by noticing that 
scenario (c) improves the fit over the standard case (case 
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model 


X 2 without LSND 


X 2 including LSND 


(a) 


1097.6 


1104.3 


(b) 


1037.8 


1044.4 


(c) 


45.7 


45.9 


(d) 


52.5 


52.7 


(e) 


53.9 


60.7 



TABLE I: x 2 obtained for (a) pure decoherence in antineu- 
trino sector, (b) pure decoherence in both sectors, (c) mixing 
plus decoherence in the antineutrino sector only, (d) mixing 
plus decoherence in both sectors, (e) standard scenario with 
and without the LSND result for all experimental data but 
KamLAND spectral distortions. 



model 


X 2 without LSND 


X 2 including LSND 


(a) 


1135.8 


1142.3 


(b) 


1095.6 


1102.3 


(c) 


83.9 


84.2 


(d) 


90.7 


90.9 


(e) 


78.6 


85.4 



TABLE II: x 2 obtained for (a) pure decoherence in antineu- 
trino sector, (b) pure decoherence in both sectors, (c) mixing 
plus decoherence in the antineutrino sector only, (d) mixing 
plus decoherence in both sectors, (e) standard scenario with 
and without the LSND result, for all data. 

(e)) by A\ 2 = 14.8 (see Table 1) but the price to pay for 
it is to give up spectral distortions in KamLAND. In fact, 
for KamLAND data alone we find a x 2 = 24.7 for case 
(e) and x 2 = 38.4 for an energy independent suppres- 
sion (case (c)), implying a A\ 2 = 13.7, which essentialy 
means that the net effect is to put our decoherence model 
back to the standad case level. 

It is important to stress once more that our sample 
point was not obtained through a scan over all the pa- 
rameter space, but by an educated guess, and therefore 
plenty of room is left for improvements. On the other 
hand, for the mixing-only/no-decoherence scenario, we 
have taken the best fit values of the state of the art 
analysis and therefore no significant improvements are 
expected. At this point a word of warning is in order: 
although superficially it seems that scenario (d), deco- 
herence plus mixing in both sectors, provides an equally 
good fit, one should remember that including decoher- 
ence effects in the neutrino sector can have undesirable 
effects in solar neutrinos, especially due to the fact that 
decoherence effects are weighted by the distance traveled 
by the neutrino, something that may lead to seizable (not 
observed!) effects in the solar case. 

Another point to stress is that, as can be easily seen 
in the figure, decoherence plays no role for atmospheric 
neutrino energies and baselines, where the effect is given 



almost totally by the oscillations driven by the mixing. 
This was somehow expected, since decoherence is known 
not to be the leading force behind the atmospheric neu- 
trinos. In fact we are taking advantage of the fact that 
for our particular set of parameters, no significant effects 
are expected in the (2,3) channels. 

One might wonder then, whether decohering effects, 
which affect the antineutrino sector sufficiently to ac- 
count for the LSND result, have any impact on the solar- 
neutrino related parameters, measured th roug h antincu- 
trinos in the KamLAND experiment |U l35|| In order 
to answer this question, it will be sufficient to calculate 
the electron survival probability for KamLAND in our 
model, which turns out to be Pp a ^>p„ Upland— -57, in 
perfect agreement with observations. As is well known, 
KamLAND is sensitive to a bunch of different reactors 
with distances spanning from 80 to 800 km. However, the 
bulk of the signal comes from just two of those, whose 
distances are 160 and 179 km. These parameters have 
been used to compute the survival probability. Never- 
theless, decohering effects, at least of the simplified type 
discussed in this work, are not able to account for any 
spectral distortion; they rather provide an overall sup- 
pression. Thus, the spectral distortion that is apparently 
being observed by KamLAND [23], tends to favor the 
standard three-generation scenario (without KamLAND 
spectral information but including LSND the \ 2 f° r mod- 
els (c) and (e) would be 45.9 and 60.7, respectively). 
Therefore, if KamLAND evidence gets stronger a simpli- 
fied decoherence model of the class discussed here would 
be ruled-out. However, this should by no means be re- 
garded as implying that decoherence models in general 
will not survive such a case. On the contrary, it is our 
belief that CPT Violating decoherence models with com- 
plex energy dependences of the decoherence parameters 
as the one discussed here, but probably of more com- 
plicated, even non linear form, stand a good chance of 
describing nature. 

It is also interesting to notice that in our model, the 
LSND effect is not given by the phase inside the oscilla- 
tion term ( which is proportional to the mass difference) 
but rather by the decoherence factor multiplying the os- 
cillation term. Therefore the tension between LSND and 
KARMEN 36] data is naturally eliminated, because the 
difference in length leads to an exponential suppression. 
Thus, while we predict a 0.24 % anti-electron neutrino 
appearance probability for LSND, the corresponding one 
for KARMEN gets only to 0.14%. And although this 
number is already below their experimental sensitivity, 
one should notice that KARMEN combines neutrino and 
antineutrino channels in their analysis, so that the ac- 
tual appearance probability is even smaller, as our effect 
shows up only in the antineutrino sector. 

Another potential source of concern for the present 
model of decoherence might be accelerator neutrino ex- 
periments, which involve high energies and long base- 
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lines, and where the decoherence L ■ E scaling can po- 
tentially be probed. This, however, is not the case. Ac- 
celerator experiments typically join their neutrino and 
antincutrino data, with the antineutrino statistics being 
always smaller than the neutrino one. This fact, together 
with the smaller antincutrino cross section, renders our 
potential signal consistent with the background contam- 
ination. Even more, in order to constrain decoherence 
effects of the kind we are proposing here through acceler- 
ator experiments, excellent control and knowledge of the 
beam background are mandatory. The new KTeV data 
[37| on kaon decay branching ratios, for example, will 
change the v e background enough to make any conclu- 
sion on the viability of decoherence models useless. After 
all, the predicted signal in our decoherence scenario will 
be at the level of the electron neutrino contamination, 
and therefore one would need to disentangle one from 
the other. 

Having said that, it is clear that althougth this very 
simplified decoherence model (once neutrino mixing is 
taken into account) does not account for all the observa- 
tions including the LSND result better than the standard 
three-neutrino model, it certainly suggests that less sim- 
plified models are worth exploring. This scenario , which 
makes dramatic predictions for the upcoming neutrino 
experiments, expresses a strong observable form of CPT 
violation in the laboratory, and, in this sense, our fit does 
not seem to give (yet) a conclusive answer to the question 
asked in the introduction as to whether the weak form of 
CPT invariance (JJJ is violated in Nature, but it certainly 
encourages further studies. 

This CPT violating pattern, with equal mass spec- 
tra for neutrinos and antineutrinos, will have dramatic 
signatures in future neutrino oscillation experiments. 
The most striking consequence will be seen in Mini- 
B00NE (2^, According to our picture, MiniBooNE will 
be able to confirm LSND only when running in the an- 
tineutrino mode and not in the neutrino one, as decoher- 
ence effects live only in the former. Smaller but exper- 
imentally accessible signatures will be seen also in MI- 
NOS |2£j, by comparing conjugated channels (most no- 
ticeably, the muon survival probability). 

We next remark that fits with decoherence parameters 
with energy dependences of the form (|ll|l imply that the 
exponential factors e Afc * in (jjjj due to decoherence will 
modify the amplitudes of the oscillatory terms due to 
mass differences, and while one term depends on L/E 
the other one is driven by L ■ E, where we have set t = L, 
with L the oscillation length (we are working with natural 
units where c = 1). 

The order of the coefficients of these quantities, 7^ ~ 
10~ 18 ,10~ 24 (GeV) 2 , found in our sample point, im- 
plies that for energies of a few GeV, which are typical 
of the pertinent experiments, such values are not far 
from 7? ~ Amfj. If our conclusions survive the next 
round of experiments, and therefore if MiniBOONE ex- 



periment |38( confirms previous LSND claims, then this 
may be a significant result. One would be tempted to 
conclude that if the above estimate holds, this would 
probably mean that the neutrino mass differences might 
be due to quantum gravity decoherence. Theoretically 
it is still unknown how the neutrinos acquire a mass, 
or what kind of mass (Majorana or Dirac) they possess. 
There are scenaria in which the mass of neutrino may be 
due to some peculiar backgrounds of string theory for in- 
stance. If the above model turns out to be right we might 
then have, for the first time in low energy physics, an in- 
dication of a direct detection of a quantum gravity effect, 
which disguised itself as an induced decohering neutrino 
mass difference. Notice that in our sample point only an- 
tineutrinos have non-trivial decoherence parameters 77 , 
for i = 1 and 3, while the corresponding quantities in 
the neutrino sector vanish. This implies that there is a 
single cause for mass differences, the decoherence in an- 
tineutrino sector, which is compatible with common mass 
differences in both sectors. If it turns out to be true, this 
would be truly amazing. 

SPECULATIONS ON THEORETICAL MODELS 
OF DECOHERENCE 

Moreover, if the neutrino masses are actually related 
to decoherence as a result of quantum gravity, this may 
have far reaching consequences for our understanding of 
the Early stages of our Universe, and even the issue of 
Dark Energy that came up recently as a result of as- 
trophysical observations on a current acceleration of the 
Universe from either distant supernovae data |4jj or mea- 
surements on Cosmic Microwave Background tempera- 
ture fluctuations from the WMAP satellite 0] . Indeed, 
as discussed in 0, H3] , decoherence implies an absence of 
a well-defined scattering S-matrix, which in turn would 
imply CPT violation in the strong form, according to the 
theorem of 14]. A positive cosmological constant A > 
will also lead to an ill definition of an S-matrix, precisely 
due to the existence, in such a case, of an asymptotic- 
future de Sitter (inflationary) phase of the universe, with 
Hubble parameter ~ y/A, implying the existence of a cos- 
mic (Hubble) horizon. This in turn will prevent a proper 
definition of pure asymptotic states. 

We would like to point out at this stage that the 
claimed value of the dark energy density component of 
the (four-dimensional) Universe today, A ~ 10~ 122 Mp, 
with Mp ~ 10 19 GeV (the Planck mass scale), can actu- 
ally be accounted for (in an amusing coincidence?) by 
the scale of the neutrino mass differences used in or- 
der to explain the oscillation experiments. Indeed, A ~ 
[{Am 2 ) 2 /M P ]M P ~ 1Q- 122 M P for Am 2 ~ lO" 5 eV 2 , 
the order of magnitude of the solar neutrino mass differ- 
ence assumed in oscillation experiments (which is the one 
that encompasses the decoherence effects, as can be seen 
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by eq.©)- The quantum decoherence origin of this mass 
then would be in perfect agreement with the decoherence 
properties of the cosmological constant vacuum, men- 
tioned previously. We note at this stage that it is possible 
to justify theoretically the above form for A cx (Am 2 ) 2 , if 
one computes the vacuum expectation value of the stress 
tensor of the neutrino field in a Robertson- Walker space- 
time background, not with respect to the mass-eigenstate 
vacuum, but with respect to a "flavour vacuum" , which 
may not be unitarily equivalent to the former. Details 
are given in |42T ] and references therein, where we refer 
the interested reader. We only note here that the issue is 
related to the old problem of mixing in a quantum field 
theory, which has been studied in detail in j^. Need- 
less to say that we by no means consider such a problem 
as completely understood to date, but we would simply 
like to point out here, and in |42|. yet another interest- 
ing feature, namely a possible link of the neutrino mass 
difference to a non zero, non perturbative, contribution 
to the cosmological constant. 

Unfortunately, at present we have no concrete theoret- 
ical model to explain the above-described decoherence in 
the antineutrino sector, with the non universal energy de- 
pendence (|llf> . Conventional quantum-gravity decoher- 
ence scenaria assume usually a universal CPT conserv- 
ing character for the decoherence parameters between 
the particle and antiparticle sectors, and in fact they 
yield too small decoherence parameters to be detected 
experimentally. For instance, in the two-level double- 
commutator model of quantum-gravity induced decoher- 
ence of j the double commutator structure in the de- 
coherence part yields the estimate: 7 ~ (AE) 2 /Mp for 
a four-dimensional quantum-gravity decoherence model. 
Such an estimate stems from the double commutator 
of the Lindblad operators [bi,bi,p] when considered be- 
tween appropriate energy eigenstates. Here, AE is the 
energy variance, which in the case of neutrinos may 
be estimated from the Hamiltonian p + m 2 /2p to be 
AE ~ (Am 2 )/E , thereby yielding a decoherence model 
with 7 oc (Am 2 ) 2 /MpE 2 , which is universal between 
particle and antiparticle sectors, and actually too small 
to be detected in oscillation experiments. 

It is possible, however, that other models yield sig- 
nificantly larger values to match our estimates for the 
decoherence parameters above. For instance, let us con- 
sider a crude model of decoherence, based on non-critical 
string theory |2^ . according to which a modified Lind- 
blad evolution of the type © can be derived as a re- 
sult of departure from world-sheet conformal invariance. 
The latter describes the interaction of string matter with 
'fuzzy' space-time singular backgrounds, entailing 'loss of 
information'. Using for the moment some technical jar- 
gon, we mention that departure from conformal invari- 
ance is formally manifested through the appearance of a 
central charge deficit in the corresponding a model |22| . 
which from a target-space viewpoint may resemble a cos- 



mological constant A. Universes with a cosmological con- 
stant are not conformal backgrounds of string theory, as 
a result of the asymptotic future Hubble horizon. In a- 
models describing the interaction of matter with foamy 
backgrounds such A does not have to be as small as the 
actual cosmological constant of our Universe, mentioned 
above. In fact it could be due to string loop corrections, 
that is higher world-sheet topologies, and this may lead 
to violations of conformal invariance of order of the string 
scale M s . 

It has been argued in that such models of non crit- 
ical string theory lead to a Lindblad type modified evo- 
lution for reduced density matrices of four-dimensional 
string matter of the form: 

d t p = i[p,H cS ]+£ : g^[g^,p\ ■ , £ = 0(M s A) (14) 

where M s is the string scale, which, in general, may be 
different from the four-dimensional Planck mass Mp ~ 
10 19 GeV, g^ is the (fuzzy) metric background in ten di- 
mensions (the critical dimension of string theory), the : 
. . . : indicate quantum ordering of the appropriate opera- 
tors, and A denotes now the ten-dimensional dimension- 
less central-charge deficit, which, if due to string loop 
corrections, it could be of order one (in natural string 
(M,) units). One may select an antisymmetric quantum 
ordering prescription, which leads to a double commuta- 
tor structure of the decoherence terms C : [g^u, [g^ , p]] ■■ 
This will be a local contribution to matter propagation, 
due to interactions of matter with space time foam, and 
thus will be dominant as compared to global contribu- 
tions due to a non-zero Dark energy component of the 
Universe jjjj . It remains a challenge of course to explain 
how such a local decoherence will not lead to large global 
contributions to the vacuum energy. 

The decoherence terms in (|14|) are essentially propor- 
tional to the world-sheet p function for the graviton in a 
'cosmological constant' space time with metric <7 M „, i.e. 
of the form /? M „ = M^ 2 R^ ~ Ag^ 0], where is the 
Ricci tensor of the cosmological constant ten-dimensional 
non-conformal string background (anO(o' = M~ 2 ) anal- 
ysis for the world-sheet /3-functions of the pertinent a- 
model suffices for our purposes). 

If we accept this scenario for time evolution of matter 
(or antimatter), then, in our case of a non-zero string- 
loop induced cosmological constant, upon projecting on 
string theory space states \g l >, where g % are background 
fields in the target space over which the matter string 
propagates (such as g^ etc.), the double commutator 
structure yields simply 22] the square of the variance 
of the metric background as a result of interaction with 
stringy matter (Ag^) 2 . This expresses essentially the 
back reaction of matter onto space time, as a result of 
the quantum properties of the foam. 

The form of this variance depends on the details of the 
model, and it is conceivable, although we have no con- 
crete model in mind at this stage, that the variances are 
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non zero only for the antineutrino sector, and such that 
they allow to write (|14fl in the form Q , where the Lind- 
blad operator is diagonal, with various entries that could 
match our estimates above. In the case of antineutrino 
decoherence, for instance, one may associate the vari- 
ance with the uncertainty in the measurement of a typi- 
cal length in the problem, which is the oscillation length. 
In our model the latter is of order L ~ Ej Aw? in both 
particle and antiparticle sectors. Thus, an uncertainty in 
Ag^ v may be associated with an uncertainty in the co- 
variant length [i^ AL 2 ~ L 2 Ag^, from which we obtain 
the estimate: {Ag^f ~ (AL 2 /L 2 ) 2 ~ (SE/E) 2 . Here, 
SE denotes an antineutrino, say, energy uncertainty due 
to the 'fuzziness' of space time, and should not be con- 
fused with the energy variance between neutrino states 
mentioned above in the context of the model of [3(J • 

It is such uncertainties that could be energy dependent 
in a peculiar way only for the antineutrino sector, as a 
result of the strong form of CPT violation due to deco- 
herence 0, e.g. SE ~ E(E/M s ) a , for some a which 
are in general different for the different decoherence pa- 
rameters, and yield the complicated energy dependence 
argued in our work above. In fact, if the conformal invari- 
ance violations in l|14l) are large only in the antineutrino 
sector, then this may provide the natural explanation as 
to why the dark energy component of the Universe is still 
small, given that the latter will be due to matter, and 
there is a huge matter-antimatter asymmetry in the Uni- 
verse. At present, however, we have no answers to such 
issues, and our remarks above should be treated with 
caution. We find it challenging, nevertheless, to try and 
interpret theoretically, in terms of microscopic models, 
the complicated energy dependence of our decoherence 
parameters (lllll . This is currently under progress. 

In any case, given that there is no full/complete theory 
of decoherence, any guess on the order of magnitude of 
the coefficients may be proven wrong. This fact encour- 
ages phenomenological analyses like ours, whose results 
can be contrasted to experimental data and may remain 
valid independently of the theoretical speculations on the 
origin and size of the 7's. 

INSTEAD OF CONCLUSIONS 



differ between neutrinos and antineutrinos. The agree- 
ment of this fit with the data was found to be better than 
in the two-generation decoherence case [i^ , where there 
was no evidence for decoherence. This conclusion de- 
pends crucially upon the inclusion of the LSND results. 
Our analysis and conclusions are therefore in disagree- 
ment with pessimistic claims made in the recent litera- 
ture |2|]|, |22| on the non-necessity of studies of quantum- 
gravity-induced decoherence in three generation models. 

It is imperative that the LSND results are confirmed by 
future experiments, such as MiniBOONE j^, because as 
we have seen above, it is indeed possible that the study 
of neutrino physics constitutes a true and rare window 
to Planck scale physics, with far reaching consequences 
ranging from cosmological considerations to the origin 
of the dark energy component of the Universe. In the 
present work, motivated by our sample point, we put 
forward a conjecture that quantum-gravity decoherence 
may be responsible for the neutrino mass differences and 
the non-zero value of the vacuum energy density of the 
Universe today. The scaling of the decoherence parame- 
ter with the energy is a crucial test of such conjectures, 
and would be a useful theoretical guide for model build- 
ing, should the results of LSND and the present work 
survive future data. 

Before closing we must stress the fact that the dynami- 
cal semi-groups approach to decoherence, examined here, 
relies on the approximation that the quantum gravity ef- 
fects imply a linear decoherence, i.e. a linear dependence 
on the reduced density matrix of matter. This may not 
be true in a complete theory of quantum gravity, and 
hence, even if we get agreement with experimental data, 
this agreement should be interpreted with great care. 
Clearly we are in need of a complete and detailed math- 
ematical model of quantum gravity, because generic and 
model independent analyses may be misleading. Never- 
theless, the possibility of verifying experimentally devi- 
ations from conventional field theoretical treatments in 
the near future, e.g. in the context of neutrino physics, 
is by itself exciting even if it cannot give a complete an- 
swer to such fundamental questions as the Planck scale 
structure of quantum space time. 



In this work we have presented a set of decoherence 
parameters which can account for all available neutrino 
data, including LSND results. By invoking a three- 
generation decoherence scenario, combined with conven- 
tional oscillations, and flavor violations, we have man- 
aged to fit the available data and to obtain decoher- 
ence parameters which could find a natural explanation 
in some theoretical models of quantum gravity. We ex- 
plained the LSND results not by assuming CPT violat- 
ing mass differences, but by invoking CPT violation in 
its strong form, namely that the decoherence parameters 
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